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Abstract: We consider the elliptic three body problem as a perturbation of the 
circular problem. We show that for sufficiently small eccentricities of the elliptic 
problem, and for energies sufficiently close to the energy of the libration point 
L2, a Cantor set of Lapunov orbits survives the perturbation. The orbits are 
perturbed to quasi-periodic invariant tori. We show that for a certain family of 
masses of the primaries, for such tori we have transversal intersections of stable 
and unstable manifolds, which lead to chaotic dynamics involving diffusion over 
a short range of energy levels. Some parts of our argument are nonrigorous, but 
are strongly backed by numerical computations. 



1. Introduction 

In the planar restricted circular three body problem (PRC3BP) two large masses 
fj, and 1 — /-t rotate on planar circular Keplerian orbits. The problem deals with 
the motion of a third massless particle, which moves on the same plane as the 
two larger bodies under their gravitational pull. This problem was considered 
by Llibrc, Martinez and Simo in [f G] for energies of solutions close to the energy 
of the libration point L^. There it has been shown that there exists a family of 
parameters {fJ-k}^?^ f° r wn i cn we have a homoclinic orbit to the libration point 
L2. What is more, it has been shown that for fi close to any of the values 
for a Lapunov orbit around L2 with energy sufficiently close to the energy of 
L2, its stable and unstable manifolds intersect transversally. This dynamics is 
restricted to a constant energy manifold and leads to a homoclinic tangle. Later 
a similar problem has been numerically investigated by Koon, Lo, Marsden and 
Ross [13], where smaller energies were considered. In such a case the chaotic 
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dynamics is extended to include homo-clinic and hctcroclinic tangles along stable 
and unstable manifolds of Lapunov orbits around both the libration point L\ 
and L2. This has been later proven by Wilczak and Zgliczynski using a method 
of covering relations and rigorous computer assisted computations in [23], for 
the case of the Jupiter-Sun system and the energy of the comet Oterma. 

All of the above mentioned results have a common feature: since the problem 
follows from an autonomous Hamiltonian, the transversality of the intersections 
and the chaotic dynamics of the system are always restricted to a constant energy 
manifold. In this paper we are going to consider the planar restricted elliptic 
three body problem (PRE3BP), where the equations are no longer autonomous 
which, means that a change of energy of solutions is possible. We will consider the 
circular problem considered by Llibre, Martinez and Simo in [ ] and generalize 
it to allow the orbits of the two larger masses /x and 1 — fi to be elliptic with small 
eccentricities e. We will treat this as a perturbation of the circular case. We will 
show that most of the Lapunov orbits around Li persist under such perturbation. 
What is more, we will show that the rich dynamics associated with these orbits 
and obtained in [16] also survives. On top of this, it will turn out that we also 
have chaotic diffusion along the energy level. In effect the dynamics of the elliptic 
problem is by one dimension richer than the dynamics of the circular problem, 
where all solutions are restricted to a constant energy manifold. To be more 
precise, we will prove the following Theorem 

Theorem 1. For a sequence of masses {^k}kLi an d f or energies close to the 
energy of the libration point Li, for sufficiently small eccentricities e of the el- 
liptic problem, most of the Lapunov orbits around Li survive and are perturbed to 
invariant tori. What is more, there exist a homoclinic and a heteroclinic tangle 
between the surviving tori which involves diffusion in energy (Such a hetero- 
clinic tangle between invariant tori is the mechanism of the so called transition 
chains). The diffusion occurs between any two of the surviving tori on an interval 
of energies of order e 1 / 2 . 

We would like to stress here, that some parts of our proof of the above theorem 
are nonrigorous, but are strongly backed by numerical computations. We will 
explain this later in this introduction after an overview of the proof. 

The diffusion between energies follows from a mechanism similar to the 1964 
Arnold's example [2]. Arnold conjectured that this phenomenon appears in the 
three body problem. The result of this paper is a small step towards a proof 
of this conjecture, but the described dynamics does not fulfill all requirements. 
First of all, prior to perturbation we do not have a fully integrable system. 
We start with the circular problem with a setting in which we already have a 
homoclinic connection between Lapunov orbits, which in our setting play the 
role of lower dimensional normally hyperbolic invariant tori. Such systems are 
referred to as apriori unstable. Secondly, and most importantly, our diffusion is 
between energies with distance of order e 1 / 2 and not of order one. The result 
of this paper should rather be viewed as a demonstration of the fact that the 
dynamics of the PRC3BP described in [16] persists under perturbation. It is 
enriched by the occurrence of diffusion in energy, but we are unable to show 
that the diffusion is large. 

Throughout some of the so far explored examples a certain pattern can be 
observed in the methods with which diffusion is detected in apriori unstable 



Transition Tori in the Planar Restricted Elliptic Three Body Problem 



3 



systems (see for example [ I 7] for a result of Moeckel on detection of transition 
tori in the case of the planar five body problem; or the work of Dclshams, Llave 
and Seara [5] on diffusion of energy for perturbations of geodesic flows on a two 
dimensional torus; also Wiggins [20], [21] discusses this mechanism in the case of 
perturbations of completely integrablc systems) . First a normally hyperbolic in- 
variant manifold foliated by invariant tori is found. The tori are required to have 
hyperbolic stable and unstable manifolds and a transversal intersection of these 
manifolds. Secondly a perturbation of the system is considered. By perturbation 
theory ([11], [22]) of normally hyperbolic manifolds, the normally hyperbolic 
invariant manifold and its stable and unstable manifolds persist under the per- 
turbation. Next step is to show that on the perturbed invariant manifold most 
of the invariant tori survive. This under appropriate nondegencracy conditions 
is a result of the Kolmogorov Arnold Moser Theorem (KAM) [3], [12]. Using 
more recent versions of the theorem (for example [10], [11] or [25]) it can be 
shown that most of the invariant tori persist and form a Cantor set having a 
positive measure in the invariant manifold. The last step is to show that the sta- 
ble and unstable manifolds of the surviving tori intersect transvcrsally. This can 
be done by the use of a Melnikov type method along a homoclinic orbit of the 
unperturbed problem. The transversal intersections between the invariant man- 
ifolds of the perturbed tori lead to homoclinic tangles for each of the surviving 
tori. In addition to this we also have a chaotic diffusion along the Cantor set of 
homoclinic tangles between the tori. In this paper we will follow this procedure. 

When applying the method to prove the existence of transition chains for a 
given physical problem the steps of the above described procedure, which present 
the biggest obstacles are usually verification of the assumptions of the KAM 
theorem and computation of the Melnikov integral. In our case the twist property 
needed for KAM theorem follows from Moser's version of the Lapunov Theorem 
[18]. For sufficiently small \x we will prove the twist property by approximating 
the PRC3BP with the Hill's problem. For larger masses fik though we need to 
back the argument by numerical computation of the twist coefficient in Moser's 
expansion. For the Melnikov integral the result is also partly rigorous and partly 
backed up by numerics. We can prove that the Melnikov integral is convergent 
and that the Melnikov function has zeros. Showing though that these zeros are 
nontrivial needs to be backed up by numerical computation of its derivative. 
We believe though that the above mentioned computations can be performed 
using rigorous-computer-assisted approach in the spirit of [23]. Such arguments 
require careful estimates, use of topological tools, and are the subject of ongoing 
work. 

The paper is organized as follows. Section two contains preliminaries, where 
we recall the earlier results on the planar restricted three body problem of [ ] , 
and introduce basic facts about the Hill's problem and the PRE3BP. In section 
three we present Moser's version of the Lapunov theorem [18] and show how 
to apply it to obtain the twist property. In section four we show that we have 
a twist on the family of Lapunov orbits around L2. In Section five we apply 
the normally hyperbolic invariant manifold theorem together with the KAM 
Theorem to show that most of the Lapunov orbits around L 2 persist under 
perturbation from the circular problem to the elliptic problem. In Section six we 
use a Melnikov type argument to detect the transversal intersections between 
the stable and unstable manifolds of the perturbed Lapunov orbits. In section 
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seven we compute the Mclnikov integral. In section eight we gather together our 
results and prove Theorem 1. 



2. Preliminaries 



2.1. The Planar Restricted Circular Three Body Problem. In the Planar re- 
stricted circular three body problem (PRC3BP) we consider the motion of a 
small massless particle under the gravitational pull of two larger bodies of mass 
p and 1 — fi, which move around the origin on circular orbits of period 2tt on 
the same plane as the massless body. The Hamiltonian of the problem is given 

by [i] 



H(q,p,t) 



P1+P2 i - M V- 



(1) 



2 n(t) r 2 (ty 

where (p,q) = (qi, ?2>Pi>P2) are the coordinates of the massless particle and 
r\ it) and r 2 it) are the distances from the masses 1 — /i and fi respectively. After 
introducing a new coordinates system {x,y,p x ,p y ) 



x = qi cos t + qi sin t, p x 
y = — qi sin t + q 2 cos t, p y 



Pi cost + p 2 sini, 
—pi sin t + p2 cos t, 



(2) 



which rotates together with the two larger masses, the larger masses become 
motionless and one obtains [1] an autonomous Hamiltonian 



where 



H(x,y,p x ,p y ) 



i2{x,y) 



(Px + y) 2 + {Py - xf 



tt{x 7 y), 



(3) 



1-M 



y 



H 

r 2 = 



y/(x + 1 - p) 2 + y 1 



The motion of the particle is given by the equation 

x = JVH(x), 



(4) 



where x = (x, y,p x ,Py) G R 4 , J 



and id is a two dimensional identity 



id 
^ -id 
matrix. 

The movement of the flow (4) is restricted to the hypcrsurfaces determined 
by the energy level C, 

M(p,C) = {(x,y,p x ,p y ) G R 4 \H(x,y, Px , Py ) = C}. 

In the x, y coordinates this means that the movement is restricted to the so 
called Hill's region defined by 

R(ti,C) = {(x,y) eR 2 \n(x,y) > -C}. 

The shape of the Hill's region R(fj,,C) will differ with C (see Figure 1). The 
focus of our attention in this paper will be on the case when the energy C is 
equal to or slightly larger than C%. For the energy C equal to C2 we have the 
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Fig. 1. The Hill's region for various energy levels: when the energy C is smaller than C£ (a), 
when C = C£ (b) and when C = Of > C% (c). 

libration point which is of the form (— fc,0,0, —k) with k > 0. The linearized 
vector field at the point has two real and two purely imaginary eigenvalues, 
thus it follows [16] from the Lapunov theorem that for energies C larger and 
sufficiently close to there exists a family of periodic orbits 1^{C) emanating 
from the equilibrium point L^- 




Fig. 2. The unstable manifold of Lj in the x, y coordinates. 



The PRC3BP admits the following reversing symmetry 

S(x,y,Px,p y ) = {x,-y,-p x ,p y ). 

We will say that an orbit q(t) is S'-symmetric when 

S(q(t))=q(-t). (5) 

We have the following results about the stable and unstable manifolds of L% 
and lfj,(C). 

Theorem 2 ([16, Theorem A]). For fj, sufficiently small the branch of 
contained in the S region (see Figures 1 and 2) has a projection on the bounded 
component of C) given by 

d(t) = /z 1/3 (^7V(oo) - 3 1/6 + M (oo) cost + o(l)), (6) 
a(t) = -tt + fi 1/3 (N(oo)t + 2M(oo)sint + o(l)), (7) 
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where d is the distance to the z.v.c., a the angular coordinate, N(po) and M(oo) 
are constants and the expressions remain true out of a given neighborhood of L^. 
The parameter t means the physical time from a suitable origin. The terms o(l) 
tend to zero when \x does and they are uniform in t for t = 0(/.i -1 / 3 ). 

In particular the first intersection with the x axis is orthogonal to that axis, 
giving a S -symmetric homoclinic orbit for a sequence of values fi which has the 
following asymptotical expression: 

Let us now introduce a notation for the S'-symmetric homoclinic to L2 orbit 
obtained in Theorem 2 for the parameters /ik given in (8). We will denote such 
an orbit by q°(t) = q1{t) (see Figure 2). Note that such an orbit starts at the 
section {y = 0} at time t = 0. 

Theorem 3 ([ , Theorem B]). For fi and AC = C — sufficiently small, 
the branch of W u (1^ (C)) contained in the S region intersects the plane y = 
for x > in a curve diffeomorphic to a circle (see Figure 3). 

What is more for points in the {pi, C) plane such that there exists a pk of 
Theorem 2 for which 

Z\C>L/4 /3 (/i-/x fe ) 2 

holds (where L is a constant), there exist S -symmetric transversal homoclinic 
orbits. In particular, for fi — fik there exist symmetrical transversal homoclinic 
orbits q^ k c f or ^ e periodic orbit lp k (C) for every C € (C2 k , Sfj, k )- 

Remark 1. In the original version of the Theorems 2 and 3 the C was taken as 
the Jacobi constant C — F where 

F(x,y,p x ,p y ) = -2H(x,y,p x ,p y ) = 2Q(x,y) - (x 2 + y 2 ) . (9) 

In this paper we have rewritten the Theorems with C as the Hamiltonian of the 
PRC3BP, which means that we have a change of sign in C compared with the 
original version. 

Remark 2. Using a standard dynamical system theory argument, from the Birkhoff- 
Smalc homoclinic theorem, the transversal homoclinic connections between La- 
punov orbits imply chaotic symbolic dynamics of the system. Since the system 
is autonomous this dynamics is restricted to the constant energy manifold. 



2.2. The Hill's Problem. Let us consider a change of coordinates which shifts 
the origin to the smaller body of the mass p, and rescales the coordinates by the 
factor /1 -1 / 3 . We will refer to the following as the Hill's coordinates. 



X 


= M" 


V3 ( x + 1 _ 


fi 


= M" 


1/3 y 


Px 


= M~ 




Py 


= M~ 


1/3 (.Py + l- 
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Fig. 3. The intersections of the stable and unstable manifolds of Z M (C) in the PRC3BP. 



Wc will rewrite the Hamiltonian (3) and derive the formula of the Hill's problem 
and list a number of facts which will be relevant for us in the future. 
Let us start with a simple lemma. 

Lemma 1. Consider a Hamiltonian H : W l x R™ — > R and a transformation 
p = \ip, q = fiq. Then (p(t),q(t)) is a solution of the Hamiltonian system for 
H(p,q) if and only if (p(t),q(t)) is a solution of the Hamiltonian system for 

The shift of the origin present in the transformation (10) is clearly canoni- 
cal hence we can use Lemma 1 to quickly obtain the Hamiltonian in the new 
variables. It is easy to see that 

(Px + y) 2 + iPy - xf = 2/3 (p x + V? + (Py - xf 
2 M 2 

By expanding Q in the new coordinates around zero one obtains the following 
family of Hamiltonians (parameterized by fi) 

H{», x, y,p- x ,p- y ) = & + y) 2 + (Py-*) 2 1 _ 3- 2 + 0(M i/3 )0(f2) + 

where C(fi) = /i -2 / 3 ^(1 - /i) + (1 — /i) 2 /2J and f = yjx 2 + y 2 . We can drop 

the term C{fi). It should be stressed that H depends analytically on x,y,p x ,p y 
and £t 1//3 , hence the derivatives of the 0(fj}^ 3 ) with respect to x, y,p x ,Py are still 
0(^ 1/3 ). Hence 

x^ 1 = JVH(n, x^ 1 ) = JVif(0, x^ 1 ) + 0(/ /3 ), 

where x = {x,y,p x ,py) and the superscript /i indicates that wc consider a Hamil- 
tonian system with fixed [i. The term O^ 1 / 3 ) is uniform for |x| < R. The 
Hamiltonian if(0,x M ) is the Hamiltonian of the Hill's problem 

ff Hill (x) = F(0,x). (11) 
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Let us now list a few of the properties of the Hill's problem. The prob- 
lem has two equilibrium points, Lf m = (-3~ 1/3 , 0, 0, -3~ 1/3 ) and Lf in = 
(S-^AO^- 1 / 3 ). The linearization of x' = JVH HlU at L% ai is given by 
x' = Ax, where 




A = 8 1 • W 
V -4-1 0/ 

The eigenvalues of A arc: ±a±, ±a2 with a\ = \/\ + 2\[l and «2 = V 7 ! - 2V7. 



2.3. The Planar Restricted Elliptic Three Body Problem. The planar restricted 
elliptic three body problem (PRE3BP) differs from the PRC3BP by the fact that 
the two larger bodies move on elliptic orbits of eccentricities e instead of circular 
orbits. The period of these orbits is 2n and the Hamiltonian of the PRE3BP is 
analogous to (1), with the only difference that in r\(t) and ^(i) we take the dis- 
tance from the elliptic instead of the circular orbits of the two larger masses. The 
trajectories of these orbits can be written as (see [24]) ((1 — /i) x (t) , (1 — fi) y (t)) 
for the body fj, and ([ix (t) , fiy (t)) for (1 — fi), where 

x(t) = (1 — ecosi/j) cosip + 0(e 2 ), 

y(i) = (l-ecosi/>)sin^ + 0(e 2 ), (13) 
ip(t) =t + 2esinf + 0(e 2 ). 

If one changes into the rotating coordinates (2) x = (x, y^VxiVy) then the Hamil- 
tonian (1) becomes 

H e (x,f) = H(x) +eG(x,t) + 0(e 2 ), (14) 

where H is the Hamiltonian of the PRC3BP (3), G is 2ir periodic over t and is 
given by the formula 

G = —— ^i(x, y, fi, t) + —^-3%, y, M - 1, t), (15) 
M (r 2 ) 

where 

rj = (x - fi) 2 + y 2 
rl = {x+l- A t) 2 + y 2 
f(x, y, a, t) = — ya[3 sini — sin 3 t] + xa[cos t + cos 3 t] — a 2 cos(i). 

The Hamiltonian H e generates a differential equation 

x = /(x)+ eff (x,t) + 0(e 2 ) (16) 

where 



/(x) = JVff(x) 
<?(x,i) = JVG(x,i). 



(17) 
(18) 
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Remark 3. The bounds 0(e 2 ) in the formulas (14) and (16) are uniform on any 
given bounded set which is separated from the two larger masses. In the future 
discussion we will be interested in the solutions of the PRE3BP which are close 
to the Libration point , to the Lapunov orbits (C) and to the orbit q^ k c 
(from the Theorem 3) for C slightly larger than C%. All of the above mentioned 
clearly satisfy this condition (see Figure 3). 

3. From the Lapunov-Moser Theorem to the Twist Property at 
Equilibrium Points 

In this section we will show how one can prove the twist property at an equilib- 
rium point using the Lapunov-Moser Theorem [18]. First the Theorem will be 
stated. Next a number of observations on the Theorem in the special case of one 
real and one pure imaginary eigenvalue with just two degrees of freedom will be 
made. This will be followed by a brief outline of the construction by which the 
Theorem was proved [18] from which the twist property will follow. 

Theorem 4 (The Lapunov-Moser Theorem [ 8] ) . Let 

x v = H Vv (x,y) (19) 
Vv = -H x „(x,y) 

v = 1, . . . ,n, be an analytic Hamiltonian system with n degrees of freedom and an 
equilibrium solution x = y = 0. Let ct\, . . . , a„, —oti, ■ ■ ■ , —a n be the eigenvalues 
of the linearization of (19) at the equilibrium point x = 0. Assume that the 
eigenvalues 

OL\ , . . . , OL n j (X \ j . . . , OL n 

are 2n different complex numbers and that a±, oli are independent over the reals. 
Let us also assume that for any integer numbers n\ and ri2 

otv 7^ niai + niui for v > 3. 

Then there exists a four parameter family of solutions of (19) of the form 

xv = <l>v{€i,b,'m,m) ( 2 °) 
y v = 4>,y(£i,&,m,V2) 

where 

Ut) = 4V a ^W«^°\ Vk (t) = V ° e -^>UW 2 ) for k = i i2 , (21) 

and 

ai&vl&ri) = «! + ..., a 2 (^?,£ 2 % ) = a 2 + ... (22) 

are convergent power series. The series <f>„,i()v converge in the neighborhood of 
the origin and the rank of the matrix 

( <?W* 4>vrt k \ 

fc=l,2,...,n 

is four. The solutions (20) depend on four complex parameters rf^. 
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In the case of the PRC3BP, n is simply equal to two and the equations (20), 
(21) describe all the solutions near the neighborhood of the equilibrium point. 
We will be interested in the application of the above Theorem to the libration 
point L2, where ct\ is real and 0,2 is pure imaginary. From now on we will restrict 
our discussion to this particular case. The following remarks and lemmas adopt 
Theorem 4 to this setting. 

Remark 4. When the system (19) is generated by a real Hamiltonian and if a\ 
is real and ot2 is pure imaginary then for the real solutions of (19) of the form 
(20) the functions and r]k(t) arc invariant under the involution [15, page 

102] 

JwiZufaviiTfr) = (£,iJm,m,i^2)- (23) 

Let us also note that the original version of Theorem 4 in [IS] contained an error. 
There an involution 

Jw{S,i^2,m,m) = (fi,»72,»7i,6) 

was proposed. This stands in conflict with a requirement that the transformation 
& in the proof of Theorem 4 [18] should be canonical (See also equations (38) 
and (39)). 

The reality condition (the fact that a point (£i, £2, ?7i, V2) given in new coor- 
dinates represents a point from IR 4 in original ones) is 

J w (£,1, £,2, 11,12) = (tufanunz)- ( 24 ) 

It is easy to see, that the subspace of C 4 of fixed points of J u , Fix(J w ) is given 

by 

6 € K, m G R. £2 = re 1 *, r, 2 = ire"'* 
where r, ip £ K. On Fix(J u) ) we will use the coordinates (£1, 771, r, ip). 

Remark 5. From the proof of the convergence of the series (20), (22) during the 
proof of Theorem 4 in [18], it follows that if we consider a family of analytic 
Hamiltonians 

H x : E" xl"^R, 

which depends analytically on a parameter A € E, then the radius of convergence 
of the series (20), (22) can be chosen uniformly for close values of A. 

Lemma 2. If ol\ is real and a.2 is pure imaginary then for all real solutions of 
(19) the series a\ from the Theorem 4 is real and the series 02 is pure imaginary. 
What is more if we choose a real periodic solution 

a?„(t)=&,(0,6W,0,»b(*)) 

where £2 (t) and r\2 (t) are given by ( 21 ), then there exist two real numbers r and 
4> such that 

&(*) = re ta2(0 ' lr2)+ ^ 

m {t) = i re -*»2(0.f 2 )-^. 
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Proof. From Remark 4 we know that the real solutions satisfy the reality con- 
dition (24). We therefore have 



^ e tai(Gvi^rS) = £0 e t ai (t°vH°2V%) (26) 

£O e ta 2 (e?„?,€ 2 % 2 °) = i (27) 
rfte-t^i&vUlvl) = ^Og-tarCe?,?,^^) (28) 
fa-ta^vUlvl) = i ^o e ta 2{i Wi^n° 2 ^ (29) 

if we choose t = then from the above we can see that ^ and 77° are real and 
that £° = *t?2- Using the fact that G K with (26) or (28) we can see that 

ai must be real. Using (27) or (29) and the fact that £2 = ^2 we can see that 
0,2 is pure imaginary. 

All periodic solutions have the initial conditions = ?/? = 0. If m addition 
we choose £2 °f the form £2 = re 1 '? then for the solution to be real, from (27), we 
must have £2 = * ? ?2- I n sucn case equation (21) gives us the periodic solutions as 

= ^°e ta2(0 '^^ ) = re ta2 ^ ir2 ^, (30) 

Lemma 2 shows that all periodic solutions of (19) which are real and lie close 
to the equilibrium point, are given by the equation 

l r (t) = (0, re ta2(0 ' lr2)+ ^, 0, ire- ta2(0 < lr2) -^), (31) 

when seen in the r\ v coordinates. Let us denote the set which contains these 
orbits by 

B R = {(0, re tf> , 0, ire' %e )\e € [0, 2tt), < r < R}, (32) 

where R is sufficiently small for the series 02(0, ir 2 ) to be convergent for r < R. 
Let P : K 4 — > R 4 be the time 2-7T shift along the trajectory of (19) i.e. 

P(q(t)) = q(t + 2tt), (33) 

where q(t) is a solution of (19). 

Lemma 3. If in the series a 2 from Theorem 4 i.e. 

02(6 m, 6^2) =a 2 + a 2 s^ir] 1 + 02,26^2 + ••• (34) 

we have 02,2 7^ 0, then for a sufficiently small R, the time 2n shift along the 
trajectory P restricted to the set Br is an analytic twist map i.e. 

P(r,») = (r, » + /(r)) (35) 
dr 
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Proof. In the £, r\ coordinates on Br from (31) we can see that the map P takes 
form 

P(0,re i6 ,0,ire- ie ) = (0, r e 2 ^^ +ie , 0, ire -^a 2 (o,ir 2 )-i0\ 

Keeping in mind that a 2 is pure imaginary we can see that P(r, 9) = (r,6 — 
i2ira 2 (0, ir 2 )). Since 02(0, ir 2 ) = a2 + a 2 , 2 'i'r 2 + 0(r 4 ) it is evident that if (12,2 7^ 0, 
then for sufficiently small r 

^(a 2 (0,ir 2 )) ^0. (36) 

We will now show how to determine whether for a given problem (19) we have 
0-2.2 7^ 0. For this we will quickly outline the construction of Moser [18] in order 
to obtain a formula for 02,2- The construction is performed in the following two 
steps 

c 4 ^ c 4 $ c 4 , 

(£i)6>m,»72) !->■ (xi,x 2 ,y u y 2 ) i-> (x,y,p y ,p y ), 

where the transformation <P changes the system (19) in the x, y,p y ,p y coordinates 
into a system with a simplified form 

x v oi iy x v + fv{x, y) \ "2 (37) 

Vu = -a v y v + g v (x,y) _ ' ' 

where tt\ and a 2 are the eigenvalues of the equilibrium point and / and g are 
power series starting from quadratic terms. From the simplified form (37) the 
transformation determines the series from Theorem 4. 

The transformation <& is a linear function which changes the coordinates so 
that the linear part of the equations (19) in the new coordinates becomes gen- 
erated by a diagonal matrix. What is more, the transformation <P should be 
canonical i.e. 

<P T J<P = J (38) 

where 

J= (-j/o) and Id= (ll 

On top of that <P should satisfy the following reality condition [15,18], which 
expresses that fact that J w is simply the map describing how the complex con- 
jugation works in new coordinates 

J Z $ = $J W , (39) 

where 

Jz(x,y,Px,P y ) = (x,y,p x ,p y ) (40) 

J w (xi,x 2 ,yi,y 2 ) = (xi,iy~2,yi,ix 2 ). 

The construction of the transformation is done by comparison of coeffi- 
cients. We look for 



v = (MS, v),<h(Z, v), MZ, v), MS, »?)), 
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with power series <f> v , ip v , a v , v = 1, 2 of the form 

<M£i;6,?7i,?72) = ELi <Wfc + h - 0± - 
ipv{€i,&,m,m) = J2l=i^kVk + h.o.t. 



(41) 



such that 



satisfy (37) if 



x v = 0„(£i,&,f/i>^) (42) 



ik = a k (Z 1 r h ,&r)2)£k (43) 
To construct (f) u ,ip u , a v one can rewrite using (42) and (43) the equation (37) as 



= eLi (iff - f% afe7 ? fe ) = °- v ^ v + /"(^ ^) 

Vu = ELl {j^ a kik - 7^ a kVkJ = -ol v i\} v + gr„(0, V), 



1,2. (44) 



and compare the coefficients in (44). Let us denote by <j) v ,N> Vv,jv, <Vjv the coef- 
ficients in the series </v, i/v, a v which come from the homogenous polynomials of 
order N. We can rewrite the part of (44) which contains all the terms of order 
N as 



Efe=l a k [£kQ^ - Vk-^; ) 4>v,N + ■ ■ ■ + S v k£kO>k,N-l = Otv4>v, 



N 



Efc=i a fe (^ fe gfr ~~ ^ afr ) ^"'-^ + ' ' ' _ dvkVkak,N-i = -a v ijjv,N + 



(45) 



where the dots indicate all the terms which can be computed from (f>i/j, ipv,l, &v,l—i 
with I = 1,...,N- 1. 

The nature of equations (45) suggests that the series can be constructed by 
induction starting with the lowest terms. It turns out though that not all of the 
coefficients can be computed from (45). This is because some of the terms in (45) 
cancel each other out. If we consider a homogenous polynomial c^" 1 ^" 11 ^ 2 ?7™ 2 
of order N from <f> v ,Ni such term will cancel out in (45) if 

^ a k U~ - m^A eff 1 ^ 1 W " a^T??* 1 W = o. 



fc=l 



or), 



This can happen only if we have 

2 



^ a fe (nfc - m fc ) - a v = 0. (46) 



fe=i 



By the assumption of the Theorem 4 that for any (GRwe have tot\ + «2 ^ 0, 
we can see that (46) is true only for the terms of the form c£, v (£,iT)i) ni (^2?72)™ 2 ■ 
The value of the coefficient c corresponding to such terms is chosen from an 
appropriate normalization condition [18]. In our case though the choice of the 
normalization does not play an important role. We are interested in computation 
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of the term 0,2,2 and this can be done by induction starting with N = 1 and 
stopping at N = 3. For N = 1 all the terms are uniquely determined. For N — 2 
there arc no terms which would cancel out in (45). For N = 3 we can use the first 
equation from (45) to find 02.2- There the coefficient 02,2 stands together with 
£2 (^2^2) and the term for £ 2 (£2^2) m 4>v,3 will cancel out from the equation. 
Therefore a 2j 2 is uniquely determined since it depends only on <^,i, Vv,i, </>i,,2 
and Vv,2- Using this procedure we can obtain a direct formula for the term 02 2- 



Lemma 4. // 



/„(afi, 0:2,3/1, 2/2) = Ei,j,fc,i>i Kiki x \ x 2yXy\ 
gu[xi,x2,yi,y2) = Eu.w^s^i^yi^ 



1,2. 



a 2,2 = ( — /l, 1,0, o/o.l, 0,1 — fo, 1.1, Offo, 1.0,1 + /o, 0,1, 150,2, 0,0 — /o,2,0,o/o. 1,0.1 (^^) 
+ 2ffo, 2, 0,o/o, 0,0, 2 + /l, 0,0,l/o, 2,0,0 — So, 1,0, l/o, l,0,l) + /o, 2,0,1 

Proof. The above can be checked from the formula (45) by direct computation. 
This ends our construction of the coefficient 0,2.2- 

Let us now briefly turn to the relation between the energy and the radius r 
of the periodic orbits (31). 

Lemma 5. For sufficiently small r the energy of the orbit l r i.e. 

h(r)=H($W r ))), (48) 

is equal to 

h{r) = H(0) + ^D 2 H(0) (<2>(0, 1, 0, »)) r 2 + o(r 2 ). (49) 

Proof. Since the problem (19) is autonomous the energy is constant along the or- 
bit l r . Without any loss of generality we can therefore assume that </> in equation 
(31) for l r is zero and compute 

h{r)=H($m r (0)))). 

Let us first note that the construction of = (<^i, 02, V'liV^) produced power 
series of the form (41), hence 

#(i r (0)) = (^^Vi^XUO)) = (0,r,0,ir) + O(r 2 ). 

The transformation <P is linear and therefore 

<P (#(J P (0))) = r${0, 1, 0, i) + 0(r 2 ). (50) 

We can compute h(r) as 

h(r)=H^m r (0)))) (51) 
= H(0)+DH(0) ($(V(l r (0)))) 

+ \d 2 h{q) (* mm)) + 0(1* (^(o)))i 2 ). 

Since zero is an equilibrium point we know that DH(0) = 0, thus by substituting 
(50) into (51) we obtain our claim. 
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4. Twist in the PRC3BP at L 2 

In this section we will show that in the PRC3BP the time 2-k shift along the 
trajectory is a twist map on the set of the Lapunov orbits around L%. This will 
be done by first proving the twist property in the Hill's problem and then using 
the Hill's problem as the approximation of the PRC3BP for sufficiently small /i. 

As mentioned in Section 2, we have a family of periodic Lapunov orbits (c) 
around for energies larger and sufficiently close to the energy C% of the 
libration point L%. This family of orbits corresponds to the set Br (see (32)) 
of orbits constructed in the previous section. In this section we will show that 
for sufficiently small /x we have the twist property on this family of periodic 
orbits. The main idea for the proof is to approximate the PRC3BP with the 
Hill's problem (11) where we have the explicit formulas for the libration point 
L% , the linearized vector field at iJf ^ s eigenvalues etc., which will allow us to 
compute the twist coefficient O22 • We will show that if aJpj 1 is nonzero then the 
same will hold for the twist coefficient for the PRC3BP with sufficiently small 
li. 

Lemma 6. Let P HM be the time 2n shift along the trajectory Poincare map of 
the Hill's problem (11). Then the P HlU is a twist map on the set of the Lapunov 
orbits which lie sufficiently close to L^ 111 . 

Proof. We will apply the procedure from the previous section and compute CI2 2 

for the equilibrium point = (3 -1 / 3 , 0, 0, 3 -1 / 3 ). The linear terms of (11) 

in L^ m are given by (12) with eigenvalues ±ai, ±0:2, ct\ = \fl + 2\/7 and 

a 2 = yi — 2\/7. The first of the two is real and the second is pure imaginary. 
We will choose the function ^ HlU composed of the eigenvectors of the eigenvalues 
±ai and ±a 2 

/ Ai 3 A 2 -0 \ 



Hill 



-9Ai 7 * v -89 1 \- \ 9A 2 7 7= s -i/39 7 7= x 

Ql(V7+4) ' a 2 (V7-4) z qi(V7+4) A a 2 (V7-4) 

9Ai-^ -9A 2 -^ 3 - T -09-^- 



au(\/7+4 
2 



(52) 



with Ai,Aj 6 1, |8 6 C. The above transformation <£ HlU satisfies the reality 
condition (39), and if we choose the coefficients Ai, A2, 3 as 



(53) 




then ^ Hl11 is also canonical. Computing the power series /" and g v from Lemma 
4 at L^ al and using (47) to compute the term aJpj , by rather laborious compu- 
tations (performed in Maple) one will obtain 



Hill V 7 ^ 



= — ( 102V7 - 57] « 1. 976 7. 
224 
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Since a^2 l 7^ 0, by Lemma 3 we have the twist property in the radius angle 
coordinates for all r such that < r < Rmn where i?Hiii is sufficiently small. 

One can also apply the procedure outlined in Section 3 to compute the coef- 
ficient 2 for the PRC3BP with parameters fif. ■ To do so one needs to compute 
the libration point , compute the expansion of the vector field at L^ k up 
to the order three, compute the linear change of coordinates <P = <P(nk) and 
compute «2 2 using Lemma 4. The numerical results of such computations are 
given in the below table. The values fXk chosen in the table are the numerical 
approximations of the series (8) from [16]. 



k 




0.2,2 


2/3 
Mfe 02,2 


2 


0.4253863522E-2 


74.94775503 


1.967649155 


3 


0.6752539971E-3 


255.7208206 


1.968237635 


4 


0.2192936884E-3 


541.7397907 


1.970039000 


10 


0.92907436E-5 


4466.568296 


1.973971883 


11 


0.68212830E-5 


5488.893874 


1.974219993 


12 


0.51549632E-5 


6616.434647 


1.974426964 


50 


0.582146E-7 


131563.3922 


1.976164111 


60 


0.336890E-7 


189459.7494 


1.976252225 


70 


0.212152E-7 


257884.2299 


1.976315175 


200 


0.9096E-9 


2105444.019 


1.976563023 



Table 1. The twist coefficient 02,2 for various masses /ife. 



The reason for which in the last column we compute /x fe 02,2 will become 
apparent during the proof of the following theorem. 

Theorem 5. For sufficiently small fx there exists an R(p) > sufficiently close 
to zero, such that the time 2n shift along the trajectory P^ of the PRC3BP on 
the set of Lapunov orbits -Bro) = {i(r)|0 < r < R(fJ-)} is an analytic twist map 
i.e. 

P(r,cf>) = (r,0 + /(r)) 

and 

J 7^0 foraIlre[Q,R(ji)]. 

Proof. It is easy to see that L2 (/i) depends analytically on /i 1 / 3 . First let us 
note that since the operator <P from our construction brings the derivative of the 
vector field to the Jordan form, the operator ^3 Body for the PRC3BP can be 

chosen close (depending analytically on /i 1 ^ 3 ) to the operator <£ Hl11 . The same 
can be said about the coefficients f v and g v from Lemma 4, since those come 
from the Taylor expansion up to the third order of the vector field at L 2 . What 
is more, from the proof of the Lapunov Moser Theorem 4 in [18], we know that 
the radius of convergence of the series from the Theorem 4 can be chosen to 
be independent from /1, for (i sufficiently close to zero. This means that the 
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coefficient a% 2 constructed for he PRC3BP will tend to the coefficient a^'j 1 of 
the Hill's problem. One thing to keep in mind though is that these coefficients 
are computed in the coordinates of the Hill's problem, where the equations of 
motion are given by (11). Comparing the radius rum in the Hill's problem with 
the radius r in the PRC3BP coordinates we can see that rum = \i x ^r. Since 
the coefficient 02,2 stands at the second order term in the series (34), the twist 
coefficient a 2 2 in the PRC3BP coordinates will be by (j,~ 2 / 3 larger than in the 
Hill's coordinates. This means that we have 

Urn M 2 / 3 4f Body = a™ 1 « 1. 976 7. (54) 

From Theorem 5 and the results from Table 1, it is reasonable to believe 
that we will have twist for all \Xk for k > 2. Let us point out that in the Hill's 
coordinates the radius of convergence can be chosen to be independent from \x. 
In the PRC3BP coordinates though, since rniii = ^^r, this radius will depend 
and decrease with fj, as = /j, 1 / 3 i?Hiii- This means that we will have the twist 
property only for orbits with a radius smaller than ^ 1 ' 3 i?HiU- 



5. Normal hyperbolicity, KAM theorem and the persistence of the 
Lapunov orbits 

In this Section we briefly recall some facts from the normal hyperbolicity theory 
and a version of the K.A.M. (Kolmogorov, Arnold, Moser) Theorem. We then 
apply the results to obtain the persistence result of a Cantor family of Lapunov 
orbits around L 2 for the perturbation from PRC3BP to PRE3BP. Our approach 
closely follows the method of [5] . We will therefore rewrite the theorems used in 
[5] and verify that their assumptions are satisfied in our particular setting. 
First let us recall the results concerning normal hyperbolicity. 

Definition 1 ([5, Al]). Let M be a manifold in R" and <l>t a C r , r > 1 flow on 
it. We say that a (smooth) manifold A C M - possibly with boundary - invariant 
under <P t is a-/3 normally hyperbolic when there is a bundle decomposition 

TM = TA®E S ®E U , 

invariant under the flow, and numbers C > 0, < (3 < a, such that for x € A 

ue^O \D$ t (x)v\ < Ce- at \v\ Vt > 0, (55) 
v E \D<P t {x)v\ < Ce at \v\ Vt < 0, (56) 

v€T x Ao\D$ t (x)v\<Ce l3 M\v\ Vt. (57) 

Theorem 6 ([5, A7]). Let A be a compact a-(3 normally hyperbolic manifold 
(possibly with a boundary) for the C r flow <L> tl satisfying the Definition 1. Then 
there exists a sufficiently small neighborhood U of A and a sufficiently small 
5 > such that 

1. The manifold A is C mm ^- S \ where n = a/ (3. 
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2. For any x in A, the set 

W° = {y€U: dist($t(y),$ t (x)) < Ge ( - a+s)t for t > 0} 
= {y€U: dist(<P t {y),<P t (x)) < Ce^~ s)t fort> 0} 

is a C r manifold and T X W£ = E x . 

3. The bundles E% are C"»™(r,r -<5) in x, where r = (a - /3)/f3, and 

W^ = {yeU: dist($ t (y), A) < Ce ( - Q+ ^' for t > 0} 
= {y G U : dist{<P t {y), A) < Ce { - f) - &)t for t > 0} 

is a C min ( r ' r °- <5 ) manifold. Moreover T X W\ = E s x . Finally 

wx=\J w s x . 

Moreover, we can find a p > sufficiently small and a (7 nun ( r > r ° _l5 ) diffeo- 
morphism from the bundle of balls of radius p in E S A to W A PI U. 

Remark 6. An analogous theorem can be stated for W% by considering the flow 

The following Theorem and two Remarks concern the persistence of the nor- 
mally hyperbolic manifold and its stable and unstable manifolds. 

Theorem 7 ([5, A. 14]). Let A C M (A not necessarily compact) be a-j3 nor- 
mally hyperbolic for the flow <Pt generated by the vector field X, which is uni- 
formly C r in a neighborhood U of A such that dist{M \U,A) > 0. Let ^ be 
the flow generated by another vector field Y which is C r and sufficiently C 1 
close to X. Then we can find a manifold r which is a' ' -0 hyperbolic for Y and 
C min(r, ri -S) dose tQ Ai where n = a/f3. 

The constants a' , /?' are arbitrarily close to a,0ifY is sufficiently C 1 close 
to X. 

The manifold T is the only (7 mm ( r >»"i-<5) norma lly hyperbolic manifold C° close 
to A and locally invariant under the flow of Y. 

The above Theorem is extended to give us a smooth dependence on the 
parameter by the following two remarks. 

Remark 7 (see [5, observation 1. page 390]). Assume that we have a family of 
flows <Pt,e, generated by vector fields X e which are jointly C r in all its variables 
(the base point x and the parameter e). Let A e be the normally hyperbolic 
manifold r from Theorem 7 for the flow <P t ,e- Then there exists a C min( - r ' ri -^ 
mapping F : A x L — > M, where n = a/ (3 and / C K is an interval containing 
zero, such that F(A, e) = A e and F(-, 0) is the identity. 

Remark 8 (see [5, observation 2. page 390]). For a family of flows <&t,e with the 
same assumptions as in Remark 7, there exists a (j*nm(?,ri-5) ^ ri _ ^ 
mapping R s : W S A X I — * M such that R s {W s A ,e) = W% e , R s {-,e)\ A = F(-,e), 

An analogous mapping R u also exists for W A . 
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Let us now turn to a quantitative version of the KAM Theorem used in [5]. 
Let us recall that a real number lo is called a Diophantinc number of exponent 
9 if there exists a constant C > such that 



C 

> — — 



for all p £ Z, q £ 



Theorem 8 (KAM Theorem [5, Theorem 4.8]). Let f : [0,1] x T 1 -> 
[0, 1] x T 1 be an exact symplectic C l map with I > 4. 
Assume that f = fa + 5 fx, where 

fo(I,i>) = (I,i> + A(I)), (58) 

A isC\ |#| >M, ond||/i|| c , < 1. 

Then, if 5 1 / 2 AI~ X = p is sufficiently small, for a set of Diophantine numbers 
of exponent 9 = 5/4, we can find invariant tori which are the graph of C l ~ 3 
functions u a , the motion on them is C l ~ 3 conjugate to the rotation u Q a by a, and 



\u„ — u 



°*\\c*- <n{ * 1/2 



< O (S 1 / 2 ) , and the ton cover the whole annulus except a set of 

measure smaller than O (M~ 1 5 1/ ' 2 ) . 

Moreover, if I > 6 we can find expansions 

u a = ul + 5u\ + tv, (59) 

with \\r a \\ cl - 4 < O (S 2 ) , and |K|| C! _ 4 < 0(1). 

All of the above results have been rewritten from [5] . Now we will apply them 
to the setting of the PRE3BP. We will first show that the set of the Lapunov 
orbits of the PRC3BP is normally hyperbolic. 

Let cj) e t s : R 4 -> E 4 be given by 

4>l s (x)=q(s + t), 

where q(-) is the solution for the PRE3BP with an initial condition q(s) = x. 
We will define the flow on the extended phase space : R 4 x R -> R 4 x R as 

$t(x,s) = (cl>l s (x),s + t). (60) 

Observe that the flow <fr\ is 27r periodic with respect to s variable, hence may 
be equivalently treated as a flow on M 4 x S 1 . This will later give us uniform C r 
estimates. 



Lemma 7. For a sufficiently small mass fi and for R(p) > sufficiently close 
to zero the set 

A = B R(tl) xS 1 , {l(r) x {t}\r £ [0, J?(/i)), t £ [0, 2^)} 

of Lapunov orbits of the PRC3BP (in the extended phase space) is a-[3 normally 
hyperbolic, where a > is close to the real eigenvalue ot\ at the Libration point 
L^ and (3 > can be chosen arbitrarily close to zero. 
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Proof. Consider the £i,T)i,r,8 coordinates from the previous section together 
with a time t coordinate. The £1,171 are the coordinates of the hyperbolic ex- 
pansion and r, <fi are the coordinates of the twist rotation around the libration 
point. We have 

M = {(Z u r h ,r,4>,t)\Z 1 ,r ll eR,r G [0,RQi)),<t> € [0,27r),tG [0,2tt)}. 
We can define 

£ u = {(£i,0,0,0,0)|£iGlR}, 

E s = {(0,771,0,0,0)1771 et}, 

TA = {(0, 0, r, <f>, t)\r G [0, i?( M )), <P G [0, 2n),t G [0, 2tt)}, 

then we will have TM = E U ®E S ®TA. The conditions (55) and (56) are satisfied 
with a coefficient a > close to the eigenvalue a\ at L 2 because the coordinates 
£1 and 771 are the coordinates of hyperbolic expansion and contraction. For suffi- 
ciently small fx the eigenvalue a.\ is close to = Vl + 2V7 ps 2. 508 3 of the 
Hill's problem. 

Let # be the flow in the extended phase space. From (31) for 
v = (0, 0, v r , V4, v t ) G TB R (^ 

and x € A we have 



\D$ t (x)v\ = 



Q,0,v r ,v<b +t^-(0,r)v r ,v t 
or 



therefore we see that the growth of derivative of D$ t (x) is at most linear in t. 
For any /3 > there exists a constant C > 0, such that for all v € T X A and all t 

\D$ t {x)v\<Ce m \v\. (61) 

We will now define the time 2-7T shift along a trajectory Poincare map and 
later apply the KAM Theorem to it. By Lemma 7 for e = we have an a- 
(3 normally hyperbolic invariant manifold for <P® of the form A = Sr(u) x S 1 . 
Let U be an open neighborhood of B R ^ . We will define time 2n shift along a 
trajectory Poincare map : U — > M 4 as 

We are now ready to apply Theorems 7 and 8 to obtain the following persis- 
tence result for the family of the Lapunov orbits. 

Theorem 9. For sufficiently small fi > and for some i?(/i) > sufficiently 
close to zero, there exists an eo((i) > such that for all e G [0, eo(fJ-)] the normally 
hyperbolic manifold (with a boundary; considered in the extended phase space) 
A = S_r( m ) x S 1 of the PRC3BP persists under the perturbation to PRE3BP with 
the parameter e, to a normally hyperbolic manifold (with a boundary) A e . What 
is more, for any e G [0, eo(//)] there exists a Cantor set € C [0, R(fJ-)] such that 
for any r G € there exists an invariant (two dimensional) torus 

l e (r) = {lt (r)\t €S 1 }, 

where If (r) is an invariant one dimensional torus under the map Pf . The 
family of tori l e (r) for r G £ cover A e except a set of a measure smaller than 
Oie 1 / 2 ). 
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Proof. By Lemma 7 we know that A = Br x S 1 is normally hyperbolic for the 
PRC3BP. Applying Theorem 7 and Remark 7 we obtain a family of normally 
hyperbolic manifolds A e locally invariant under and a function F : A x 
[0, e (/x)] R 4 x S 1 such that 

F(A, e)=A e = {(A^, t) \A he c M 4 , t G S 1 }. 

From the Implicit Function Theorem we know that the libration point 1)^ contin- 
ues for small values of e to a 2w periodic orbit Li2' e (t). We can modify F so that 
F(L%,t,e) = L^ e (t). By Remark 7 the function F is C ri , where n = a//3. 
From the proof of Lemma 7 we know that for sufficiently small [i we have 

a w yl + 2\/f and that /3 > can be chosen to be arbitrarily close to zero. This 
means that for sufficiently small /x, the function F is C fc for any given k > 0. 
Since yt e is locally invariant under and the flow is 2n periodic, for any to G S* 1 
the manifold vl t0je is locally invariant under P t e Q . 

Let us fix <o — and fix a Poincare map P e := Pt 0= o (here we could consider 
a map P t e for any ^ G S 1 , but we fix to = for simplicity). We will now show 
that for sufficiently small /i, e and sufficiently small < R(fJ.) the Poincare map 

P e :F(B R(M) ,0 1 e)^A Q , e (62) 

is properly defined and symplectic. For sufficiently small fi, e and sufficiently 
small R(n) we can see that the map (62) is properly defined because A 0i e is locally 
invariant. The map P e is a restriction to Ao. e of a time 2ir shift along a trajectory 
of a Hamiltonian system. Such shift is symplectic for the standard form uj = dx A 
dp x + dy A . In order to show that P e is symplectic it is therefore sufficient to 
show that uj is non degenerate on Ao, e - For sufficiently small [i and e the manifold 
ylo, e is arbitrarily close to the manifold of the Lapunov orbits of the Hill's problem 
(II), which in turn, for r sufficiently close to zero, is arbitrarily close to the vector 
space V given by the eigenvectors of the pure complex eigenvalues ±aJp n = 

±v 1 — 2\/7. To show that uj is not degenerate on Ao, e it is therefore sufficient 
to show that uj is not degenerate on V. The eigenvectors corresponding to zta^ 1 " 
are v and —iv, where v is the second column in (p HzU ; which was symplectic, 
therefore uj(v, —iv) = 1. The space V is spanned by x\ and X2, where v = xi+ix2- 
An easy computation shows that 1 = u)(v, —iv) = —2u>(xi, x%), which means that 
uj is not degenerate on V. 

Now we will use the KAM Theorem 8 to show that most of the Lapunov 
orbits on Aq iE survive under a sufficiently small perturbation. Let us first note 
that even though the Theorem is stated for a map / on [0, 1] x T 1 , the KAM 
result is local by its nature and also holds for a map / : [0, 1] x T 1 — > M x T 1 , as 
will be the case in our setting. Let uj denote the standard symplectic form in M 4 
i.e. uj = dxAdp x +dyAdpy, where (x,y,p x ,p y ) G M 4 . Let uj e denote the induced 
form on Aq, s . There exist C Tl ~ 2 (jointly with the parameter e) close to identity 
coordinate maps c e : Ao_ e — » Ao t o which transport the symplectic forms uj e into 
the standard one (see [5, page 367]). Taking r(/j.) < R((i) the map 

P e := c e oP e o (c e ) _1 : B r{fl) -> B R{fl} 

is properly defined for sufficiently small e. Clearly for e — we have P a = P . 
From the fact that P e is symplectic and the fact that P° = Pt^b is a twist map 
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follow the same properties for our maps P e and P° respectively. Now we pass 
to the action angle coordinates. From Lemma 3 we have that P° has the form 
(58). Exact simplecticity of P e is a direct consequence of simplecticity combined 
with invariance of the origin. To apply the KAM Theorem 8 what is now left is 
to show that 

\\P e -P°\\ cri ^=0(e). (63) 

This comes from the fact that in the neighborhood of the perturbing term 
e<j(x, t) + 0(e 2 ) in (16) is uniformly 0(e) in the C l norm. Observe that this esti- 
mate holds both in the original coordinates and in the action angle coordinates, 
because the origin is the fixed point for P e . Therefore the time 2tt shift maps 
Pf and P® are also O(e) close, from which (63) follows. This gives us a Cantor 
family of invariant tori (r) for r € £. Now for r £ £ we can define 

I»(r) ;={$ e t (x,0)\xel e (r),t£ [0,2n)} 
W :=^ e Go(r),0). 

We can now set R(fi) to be the radius of one of the surviving tori. This will 
ensure that A e is an invariant manifold with a boundary / e (i?(/x)). The fact that 
the complement of the Cantor set £ is 0(e 1 ^ 2 ) follows from the KAM Theorem 
(see also [19] for more details). 

Remark 9. An identical argument to the above proof can be performed to obtain 
a mirror result to Theorem 9 for any fixed parameter /Ltfc. To do so though one 
would have to verify the twist condition. This has been verified numerically for 
a sequence of parameters in Table 1. It is visible that as the masses Hk decrease 
the twist coefficient increases to infinity. It is therefore reasonable to believe that 
the twist condition holds for all parameters 

Remark 10. Let us fix a small \i and some arbitrary k > 0. Then for any e > 
there exists an interval I C [0, R((ji)] of the measure of order e 1 / 2 , for which the 
set i" n <£ for which Lapunov orbits persist under perturbation has gaps smaller 
than Ke. 

Later on this fact will allow us to construct transition chains between any two 
energies from I n £ , and is the underlying reason for which we have transitions 
of order e 1 / 2 in Theorem 1. 

Proof. The fact that such an interval exists will follow from the fact that the 
complement of the Cantor set £ is of the measure 0(e 1//2 ). Let us divide the 
interval [0, R(fJ-)] into n equal parts. If on every interval the set £ contains gaps 
larger than Ke, then from the fact that the measure of the complement of £ is 
0(e 1/2 ) (let us say that this (^(e 1 / 2 ) is equal to Me 1 / 2 for some M > 0) the 
number of such intervals n must satisfy 

une < Me 1/2 , 

which means that n < —^-pi. If we divide the interval [0, i?(/u)] into a slightly 
larger number n of equal intervals then at least one of them (this will be our 
interval I) cannot contain a gap larger than ne. The size of such an interval is 
equal to 

n ~ M 
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6. Melnikov method 

In the previous section we have shown that the normally hyperbolic manifold 
with a boundary A — B R ^ x S 1 (considered in the extended phase space) of 
the PRC3BP persists under perturbation to A e which is a normally hyperbolic 
invariant manifold with a boundary of the PRE3BP with eccentricity e. What is 
more, we have shown that A e contains a Cantor set of two dimensional invariant 
KAM tori. In this section we will consider the problem of intersections of the 
stable and unstable manifolds of such tori. 

It will be convenient for us to parametrize the manifolds A and A e using 
the radius angle coordinates r, of the Lapunov orbits from Section 3 together 
with time t G S . For e = we thus use a natural parametrization of the 
Lapunov orbits by their Birkhoff normal form coordinates (coordinates obtained 
from Theorem 4). After the perturbation it will be enough for us to use the 
fact that the parametrization is smooth (in fact, from the proof of Theorem 9 
we know that it will be C 1 * 1 with n = a/ (3) and that we can parametrize (see 
proof of Theorem 9) the perturbed libration point 1/2 (which continues to a 2tt 
periodic orbit) by r = 0. We can not assume though that the surviving perturbed 
KAM tori are parametrized by r. Our parametrization simply follows from the 
Normally Hyperbolic Invariant Manifold Theorem (see Theorem 7 and Remarks 

7, 8) without involving the KAM Theorem. 

Let us recall that prior to the perturbation the fibres W£ for p £ A intersect 
transvcrsally with the section {y = 0} (see Section 2, Theorems 2, 3 and also 
[ ]). The same goes for the unstable fibers W£. This means that for sufficiently 
small < e the same will hold for the stable fibres Wp' e and unstable fibres W"' e 
of points p G A e . Each point p G A e can be parametrized by r, <j> and t . The 
fibres of such points Wp' e , Wp' e are one dimensional and contained in sections 
S to = {(q,t )\q G K 4 }. The intersections of W^ e and of W^ e with {y = 0} can 
be smoothly parametrized by r, 4>, to and e. For a point p G A e parametrized by 
r, (j> and to we introduce the following notation for the first intersections of Wp ,e 
and of W^ e with {y = 0} 

(p s (r,cf>,to,e),t ) =w;' e n{y = 0}, 
(p u (r, ( jy,t ,e),t o )=W^ e n{y = 0}. 

Let q s (r,(f>,to,e,t) and q u (r,(p,to,e,t) be the orbits (considered in the stan- 
dard (not extended) phase space) of the PRE3BP, which start from the points 
p s (r, 0, to, e) and p u (r, 0, t , e) respectively at time t = to i.e. 

q s (r, (/), t , e, t ) = p s (r, 0, to , e) , 
q u (r, <p, t , e, t ) = p u (r, 0, t , e). 

Let us note that for parameters /x = \ik from Theorem 2 

<7 S (0, 0, t , 0, t) = <f (0, 0, t , 0, t) = q°(t - t ), (64) 

where q°(t) is the homoclinic orbit to in the PRC3BP. 
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Lemma 8. For r > and i G {s, u} we have 

9*(r, </>, to, 0, t ) = ?*(0, <£, to, 0, t ) + 0(r) 

q (r, 0, t , e, t ) = g*(r, 0, t , 0, t ) + e "^( c > 0> *o, 0, to) + o(e) (65) 

■^■(r.^to.O.to) = ^(O,0,t o ,O,t o ) + O(r), 

where the bounds o(e) and O(r) are independent from to- 
In addition 

^(^(M>*o>M)) = J D/(g i (O,0,t o ,O,t)) ^(0>,t o ,0,t) (66) 

+ .o(a*(O,0,t o ,O,t),t), 

and ^-(0, 0, to, 0, t) is bounded for all t G [to,+oo) /or i ~ s (or t G (— oo,to] 
/or i = u). 

Proof. Equations (65) come from the Taylor expansion and the fact that both 
the normally hyperbolic manifold and the foliation of its stable and unstable 
manifolds behave smoothly under the perturbation. What is more because (16) 
is 2-7T periodic in t the bounds in (65) can be chosen independent from to. 
Let i = s. For (66), using (16) we can compute 

f (O,0,t o ,O,t)) =1-1^(0,^0,0,*) 



dt \ de J de dt 

Q 

= Q- e (f (9 S (0, 0, to, e, t)) + eg (q s (0, d>, t , e, t), t) + 0(e 2 )) | e=0 

= Df (q s (0, (f>, t , 0, t)) ^1(0, 0, t , 0, t) + g (q s (0, 0, t , 0, t)) . 
oe 

For t tending to infinity the orbit g s (0, cj>, to, 0, t) converges to the point 1/2, 
and the orbit q s (0, (f>, to, e, t) converges to the periodic orbit which continues from 
1/2, hence using Gronwall type estimates we obtain 

|q s (O,0,t o , e ,t)-q s (O,0,to,O,t)| =0(e). (67) 

From Taylor expansion, for t G [to, +00) 

q s (0, 0, t , e, t) - q s (0, 0, t , 0, t) = e-^-(0, 0, to, 0, t) + re (e, t) , (68) 

where re is such that for each t wc have lim e ^ * ^ = 0. From (67) and (68) 
follows that 



^(0,<Mo,0,t) 



= - \q s {0, </>, to, e, t) - q s {0, <j>, to, 0, t) - K (e, t)| 
e 

re (e, t) 



e 

for some M > 0. Passing with e to zero gives us that ^-(r, </>, to, 0, t) is bounded 
on [to, +00). 

For i = u the argument is analogous. 
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We now have the following lemma regarding the energy of the points p s (r, 0, to , e) 
and p"(r,0,t o ,e). 

Lemma 9. Assume that is one of the parameters \Xk for which a homoclinic 
orbit q 1 (t) to exists (See Theorem 2 ). For any two points pi,P2 € A e with 
coordinates 7T,0i,io ond r<i , 02 , to respectively, we have 



ff(p s (n, X , t , e)) - H(p u (r 2 ,h, t ,e)) 

= H(l(n)) - H(l(r 2 )) + eM(to) + 0(ei?(/i)) + o(e), 



whe 



r+<x> 

M(to)= I {H,G}(q°(t-t ),t)dt. 



Proof. Let • denote the scalar product and let A s and A u denote the following 
functions 

A s (t, t ) := WH(q°(t - t )) ■ 2£(0, 0i, t , 0, t) 
A u (t, t ) := VH(q a (t - t )) ■ ^(0, 2 , t Q , 0, t). 

Using the facts that H = H r = H(l(r)) is constant along the solutions q s (r, 0, to,0, t) 
of the PRC3BP, from (65) and (64) we can compute 

H(q s (ri, 0i, to, e, t )) 

=H(q s {n, 0i, t , 0, t )) 

+ eVff (<f (r x , 0i, t , 0, to)) ■ (n , X , to, 0, t ) + o(e) 
=# ri + eVif (<z°(0) + O(n)) ■ (^(0, 0i, to, 0, t ) + 0(n)J + o(e) 

=ff ri + eVH(q°(0)) • -^-(0, 0i, to, 0, t ) + O(erj) + o(c) 
=^ ri +e^ s (t ,to)+0(eri) + o(e), (69) 

and similarly one can show that 

{q u (r 2 , 2 , t , e, t ) ) = H r2 + eA u (t , t ) + 0(er 2 ) + o(e). (70) 

Let us stress that, by Lemma 8 we know that O(eri), 0(er 2 ) and o(e) are uniform 
with respect to to- 

Let us investigate the evolution of A s (t,t ) and A u (t,t ) in time. Let us 
concentrate on the term A s (t, to)- Using (66), (64) and Vi? = —Jf (see (17)) 
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we can compute 

~ {AS, h)) = (JDf(q°(t - t ))q°(t - to)) ■ ^(0, <f>i,t Q , 0, *) 

+ (J/(g°(*-*o)))-^^(0,^,*o,0,i) (71) 
= {JDf(q°(t - t Q ))f(q°(t - to))) ■ ^(0, fc, t , 0, t) 

+ (Jf(q°(t - t ))) ■ (l>f (q°(t - t )) ^(0, 0!, t , 0, t) 

+ {Jf( q °(t-to)))-g(q (t-to),t) 
= (Jf( q °(t-to)))-g(q (t-to),t) 
= -{H,G}(q°(t-t ),t), 

where the third equality comes from the fact that for any p, q £ R 4 

{JDf(q°(t - io))p) • q + (Jp) • (Df (q°(t - t )) q) = 0, (72) 

with p = f(q°(t - t )) and q = ^(0, X , t , 0, t). Equation (72) follows from 
the fact that ui(p, q) = Jp ■ q is the standard symplectic form which is invariant 
under the flow 4>(t, x) of (4) i.e. 

( d d \ 

W I dx^ t,X ' P ' dx^ t,X ^ q ) ( 73 ) 

hence by differentiating (73) with respect to t and setting t = we obtain (72). 
We can now compute A s (to,to) using (71) 

r+oo 

A s {+w,t ) - ASo,t ) = / {H,G}(q°(t-t Q ),i)dt. 

Since limt->+oo q°(t—to) = L2 and f{Li) = 0, from the fact that ^-(0, </>i, to, 0, t) 
is bounded on [to, +00) we have 

A s (+oo,t ) = lim Jf(q°(t - to)) ~(0, <£i, t , 0,t) = 

and therefore 

/■+OC 

-ASo,to)= {H,G}(q°(t-t ),t)dt. (74) 



Analogous computations give 

r t 



ASo,t )= I {H,G}(q°(t-to),t)dt. (75) 

00 



Using the fact that r\,T2 < R(n) from (69), (70), (74) and (75) we obtain our 
claim. 



Transition Tori in the Planar Restricted Elliptic Three Body Problem 



27 



X 



e = 



X 



e < e(R) 



p I (L 2 ,f ,0) = p"(L 2 ,%0) = ? °(0) 




v J (f ,e) 




X 



X 



v"(t ,e) 



Fig. 4. The manifolds Wf, W, 1 } W.t , , and Wit, , for r > _R, intersected with XV., and 
{?/ = 0}, and projected onto the x,x coordinates. 

Theorem 10. // we consider the PRE3BP with a sufficiently small parameter 
/i = /ife for which a homoclinic orbit to exists (See Theorem 2 ) and if 



has simple zeros then for < sufficiently close to zero and any R £ 

(0,i?(^)) there exists an en(i?) such that for all e £ [0,eo(i?)] and all r £ 
£fl [i?, for which l(r) survive under perturbation and is perturbed to 



the manifolds W^,^ and Wji^ (considered in the extended phase space) inter- 
sects transver sally. 

Proof. In [] (i] it has been shown (see also Theorem 3 in Section 2.1) that for the 
unperturbed PRC3BP Wf/ r \ intersects transversally with Wu r ) at {y = 0}. Let 

U to = M 4 x {t } be the time t section in the extended phase space. Let v°(to) 
denote a point in Wf, r \ n n {y — 0} n S to ■ There can be more than just one 
such point (see Figures 3, 4), namely, if we consider Tr Xt ±(W^ n Sf n {y = 0}) 
and 7r a) i(W^ r -j fl £t H {y = 0}) then the two sets are homeomorphic to two 
circles, which intersect transversally at least one point (x = xq,x = 0) (see 
Theorem 3). Fixing any one of such points of intersection will be sufficient for 
our proof. From the construction we have Tr x ,x,y,t(v°(to)) = (xn, 0, 0, to). The 
extended phase space of the PRC3BP is five dimensional. We can choose the 
energy H to be the remaining fifth coordinate in the neighborhood of v°(to). In 
these coordinates 



Now we perturb from the PRC3BP to the PRE3BP. Let us consider some 
r, R £ M. such that < R < r < R(fi) and r £ € for which the Lapunov orbit of 
radius r is perturbed into a two dimensional (when considered in the extended 
phase space) invariant torus l e (r). Consider now the sets 




(76) 



m = m (r),to)\t e s 1 } 



v°(t ) = {x = x ,y = 0,x = 0,H = H(l(r)), t = t ). 



[w]t {r) n Eto n {y = o}) and tt x , x (wf. (r) n s to n {y = o}) 
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For sufficiently small e < e(R) these sets will remain homeomorphic to circles and 
intersect at some point (x = Xo(toi e), x = io(to, e)) which is close to (xq, 0). For 
the PRE3BP the energy is no longer preserved. This means that the intersection 
of the circles on the x, x plane does not imply an intersection in the extended 
phase space. Namely we have two points v s (t , e) <E Wf e , r \ H St H {y = 0} and 
v u (to, e) £ Wje( r j H 27t n {y = 0} which may differ on the energy coordinate (see 
Figure 5) 

v s (t„, e) = (x {t , e), 0, x (t , e), /i s (t , e), to) (77) 
u u (i , e) = (x (*o, e), 0, x (t , e), h u (t , e), i ) • 

We will show that the assumptions of our theorem imply that for some to the 
points v s (to-, e) and v u (to, e) coincide. This will imply intersection between W* e / r j 
and W£( r \. Later we will also show that such intersection is transversal. 

The points v s (to,e) and v u (to,e) are contained both in {y = 0}. What is 
more v s (t Q ,e) £ W^ (r) PI S to = W^ {r) and v u (t ,e) £ W£ (r) n = (r) . 

This means that there exist r s , (j> s and r", 0" (these depend on to and e but we 
omit this in our notations for simplicity) such that 

v s (t ,e) = (p s (r s ,(j> s ,to,e),to) (78) 
v u (t ,e) = (j, u (r u ,<t> u ,to,e),t ). 

In the proof of Theorem 9 lt a ( r ) is constructed from continuation along trajec- 
tories of a KAM torus Zg( r )- Therefore from (59) in KAM Theorem 8 we have 
that 

r s = r + 0(e), (79) 
r u =r + 0(e), 

and the bound O(e) is uniform for all r £ £. Applying Lemma 9 we have that 

h u (t , e) - h s (to, e) = H(v u (to,e)) - H(v s (to,e)) 

= H (p s (r s , s , to , e) ) - H(p u (r u , cf> u , t , e) ) (80) 
= i?(;(r s )) - H(l(r u )) + eM(to) + 0(eR(fx)) + o(e). 

The orbits ?(r s ) and i(r") are the Lapunov orbits of the PRC3BP. The radius 
of these orbits in the Birkhoff normal form coordinates is r s and r u . By Lemma 
5 we know that 

H(l(r)) = H(L 2 ) + \d 2 H(L 2 ) («P(0, 1, 0, i)) r 2 + o(r 2 ) 

hence from (79) and (80) we have 

h u (t , e) - h s (t , e) = 0(re) + eM(t ) + 0(eR(^i)) + o(e) (81) 
= eM(t ) + 0(eR(fi)) + o(e). 

Setting first R(fi) sufficiently small and then reducing e sufficiently close to 
zero implies that since M(to) has simple zeros, for some parameters to (close to 
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these zeros) we will have h u (to, e) — h s (to, e) = 0, which implies that v s (to, e) = 
v u (to, e) and in turn ensures that 

Now we will show that this intersection is transversal. First note that from 
the analyticity of the functions v s (to,e) and v u (to,e) using the same argument 
as in the proof of Lemma 9 we also have that 

- (h u (t , e) - h'(to,e)) = e^M(i) + 0(eR(ji)) + o(e). (82) 

We know that prior to our perturbation W^ r \ and W^ r \ intersect transversally 

at points v°(ta). This intersection is not transversal in the full extended phase 
space, it is only transversal in the constant energy manifold 

M = {(x,y,x,H,t)\H = H(l(r))} C R 3 x {H(l(r))} x S 1 . 

To be more precise, we know that v°(t ) £ Wf , H Wp , -,, that Wf , ■., 

W,o w C I7 t0 and that [16] 

^o (to) (Wg (r) ) +T„o (to) (WS (r) ) = M 3 x {0} x {0}. 

tQ \ / tQ V / 

These properties are preserved under small perturbation e > 0, hence for v s (to, e ) : 
v u (to, e) =: v(t Q , e) we have 

T v( t ,e)m tQ [r)) +^( i0 ,e)(^| o(r) ) = R 3 X {0} X {0}. 

We need to show that we also have transversality on the to and energy coordi- 
nate. For a fixed e the curves v s (t,e) and v u (t,e) belong to W£( r ) and W^ r ) 

respectively. At the time t = to for which i> s (io, e) = w"(<o, e) = i>(icb e) we have 
M'(t ) 7^ 0. This means that using (82), for sufficiently small e, 

^ (7T H (« s (i, e)) - n H (v u (t, e))) | t=to = | (ft«(f, e) - h s (t, e)) | t=to 

= eM'(fo) + 0(eR(fi)) + o(e) 

We also have for i e {u, s} 

_(^( t , e) ) = _ i = 1 . 

This since f t v s (t, e)\ t=to £ T v[to ^{W^ {r) ) and e)| t=to € T„ (to , e) (WJ2 (p) ) 

implies transversality, which finishes our proof. 

The order of choice of parameters in the above argument is important, so let 
us quickly run through how it should be conducted. We first choose sufficiently 
small R(fi) so that Theorem 9 and Remark 9 can be applied, and such that the 
terms 0(eR(fj)) in (81) and (82) are small. We then choose any R £ [0, R(n)] and 
choose e sufficiently small so that for any r £ [R, R(fM)] H <£ we have intersections 
of ir x ,±{W]i (r) n Sto n {y = 0}) and 7r Xli (W,'. (r) n S t0 n {y = 0}). The parameter 
e needs also to be sufficiently small so that M (to) and M'{to) dominate in (81) 
and (82) respectively. 
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h\t 2 ,e) 
h u (t 2 ,e) 



'He) 
u (t 2 ,e) 



K ( r) 



wi 
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Fig. 5. The manifolds Wjlf r % an d ^;e( r ) intersected with St 1 , St 2 an d {?/ = 0}, m the x, x, H 
coordinates for t\ < to < 



Corollary 1. J/ £/ie Melnikov integral has a simple zero then for sufficiently 
small i?(//) there exists a k > 0, such that for any two radii r± and r2 from 
£ fl [R, for which the Lapunov orbits persist under perturbation, and such 

that 

\r 1 - r 2 \ < Ke, 

the manifolds Wf e ^ r .^ and Wjf^^ intersect transversally for i,j £ {1, 2}. 

Proof. The proof of this fact is a mirror argument to the proof of Theorem 10. 
Below we restrict our attention to pointing out the difference we have connected 
with the derivation of (81) in the setting where we have two radii. 

Let /i" (in, e) and h 2 (tQ, e) stand for energies of points of potential intersection 
(constructed analogously to h u (to,e) and h s (to,e) in (77)). Let r" and r| stand 
for radii constructed analogously to r u and r s (see (78)), but coming from the 
unstable and stable manifold of l e (r\) and l e {r2) respectively. Using Lemma 5 
and the fact that r" — n = 0(e) and r 2 — r 2 = 0(e) we have 

\H(l(r?)) - H(l(r' 2 ))\ < 0(1 W) 2 - (^) 2 |) + "(%) 2 ) 

<0(|(ri) 2 -(r 2 ) 2 |) + 0(ei?( M )) + (i?( M ) 2 ) 
< 0(R(p)\ri - ra|) + 0(eR{n)) + o{R{^f). 

Using an identical argument to the derivation of (81) this gives us 

K(t ,e) - h s 2 {t ,e) = eM(t ) + O(fl0*)|r x - r 2 |) + O(efl0*)) + o(e). 

Using this estimate and following the proof of Theorem 10 we obtain our claim. 

Remark 11. Mirror arguments to the proof of Theorem 10 and Corollary 1 give 
transversal intersections of invariant manifolds for invariant tori of the PRE3BP 
with /i = fj,k for k = 2,3,.... Here we state this as a separate remark since 
Theorem 10 and Corollary 1 arc fully rigorous and do not rely on any numerical 
computations. For the argument with an arbitrary (j.^ we would need to use the 
fact that in the PRC3BP we have a twist property on the family of Lapunov or- 
bits and apply Remark 9. The twist for arbitrary fj.^ has only been demonstrated 
numerically (see Table 1). 
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7. Computation of the Melnikov integral. 

In this section we will demonstrate that for <o — and for all parameters fik from 
Theorem 2 the Melnikov integral M (t ) (76) is zero and also that ^ff(O) ^ 0. The 
fact that the Melnikov integral is zero will follow directly from the S-symmetry 
(5) of the homoclinic orbit q°(t). The fact that ^r-(O) ^ will be demonstrated 
numerically. We will compute the integral for the first few parameters fi^ and 
then demonstrate that for sufficiently small parameters fi^, the integral converges 
to an integral along an unstable manifold of the Hill's problem. 



7.1. The integral and its derivative at to = 0. We start with a lemma which 
ensures the convergence of the Melnikov integral (76). 

Lemma 10. The Melnikov integral (76) is convergent. The Melnikov function 
can be expressed as 



dt. 



id also 



dM 
~dt 



{to) 



O^G / n / \ \ O^G , _ . 

(q° (t),t + t ) - — (L 2 ,t + t ) 



dt 2 



dt 2 



dt. 



(83) 



(84) 



Proof. The orbit q° (t) is the homoclinic orbit to the Libration point L 2 . Let us 
note that the velocity x and y of q°(t) exponentially tends to zero as t tends 
to plus infinity and minus infinity. What is more the partial derivatives of G on 
q°(t) are uniformly bounded. This means that the integral over 



\{H,G}\(q°(t),t + t )dt 



, . 0G , dG , , n / \ \ -. 

1^-5- +Vs-\(9 (*)»* + *o)«ft, 



dx 



dy 



is convergent uniformly with respect to to. 

The orbit q° (t) is the solution of the PRC3BP, hence differentiating gives 



dG 
~dt 



dG 



(q°(t),t + t Q ) = ^{q° (t),t + t )+{G,H} (q°(t),t + t ) 



(85) 



From (85) we have 



M(t ) 



+ 00 



{H, G} {q°(t),t + t )dt 



= T 1 ™ (j* ^ W < 1 + **>) - § («" W ' * + ^ * 

= lim [G (g° (-T) , -T + 1 ) - G (<7° (T) , T + to) 

1 — >oo 

/ T dG 
— (q° (t),t + t )dt 



(86) 
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From the fact that liiriT^±oo q° (T) = L2 it follows that 




(87) 

uniformly with respect to to . 

From (86) and (87) we obtain (83). 

To prove (84) it is enough to observe that the formal integration of (83) is cor- 
rect, because the integral on the right hand side of (84) is uniformly convergent 
with respect to to for the same reasons as the integral in formula (83). 

It turns out that the computation of the Melnikov integral is not the ma- 
jor obstacle. The fact that we have a zero for to = follows directly from the 
^-symmetry (5) of the homoclinic orbit q°(t). This fact is shown in the below 
lemma. Later though wc will need to show that this zero is nontrivial by comput- 
ing tt(0)> which turns out to be a much harder task. The lemma also provides 
the formula for the needed integral. 

Lemma 11. The Melnikov integral (76) at to = is equal to zero and 

^'^(Sv^-^M*- (88) 

Proof. This result follows from the fact that both q° (t) and L2 are S-symmetric 
with respect to the symmetry (5). By direct computation one can check that 

$f (S (•) , -t) = -§f (•, t) and that (S (•) ,-t) = %£ (-,t) hence we have 

£.(§wm- 4 >-§H*- 

=r(^ ( '° ( -' , '-' , -^'- t) ) d< 
= ir <5 (q ° {t)) ■ _t) ■ ^ <s <i2> • _<) ) dt 

f + " (dG , „, . , dG ,, ,\ , 

= -J (at-^W^-ar^'))*' 

which gives M(0) = 0. The formula (88) follows from an analogous computation 
using 0(S(-),-t) = 0M). 

The formula for needed for the computation of (88) follows from (15) 
d 2 G 1 — u fi 

-579- = ——^ftt{x,y^,t) + — ^/ tt (ar,j/,|U- l,t), (89) 
m (n) (r 2 ) 

where 

ftt(x, y,a,t) = ya^[sint + sin3t] — xa^-[7 cost + 9cos3i] + a 2 cos(i). (90) 
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Remark 12. The verification of the fact that ^jjf(O) is nonzero is not straight- 
forward. In this paper we will restrict to numerical verification of this fact. We 
believe that for a given fx^ such computation can be performed by the use of 
rigorous-computer-assisted integration along the homoclinic orbit. Such compu- 
tations though are far from trivial and will be the subject of forthcoming work. 

In Table 2 we enclose the numerical results for the computation of (88) ob- 
tained for k up to 13. 



Table 2. Numerical results for the derivation of ^jr(O) for various mass 

parameters 

Remark 13. For the numerical computation of ^M- in (88), it is worthwhile to 
make use of the fact that the second term is 2ir periodic and that its integral 
over the period is zero. Most of the integration is done along the homoclinic orbit 
far from the region of L2, where the first term is very small and this property 
can be exploited to reduce the error of computations. The second term however 
comes in handy in a close neighborhood of L^. There it cancels out with the first 
term and therefore needs to be included. 

7.2. Approximation by the Hill's problem. The results in Table 2 seem to be 
convergent up to a sign. This is not by accident. For sufficiently small ll within a 
small region of L2, the trajectory q° can be approximated by the unstable orbit 
of L2 for the Hill's problem. This fact has been heavily exploited for the proofs of 
Theorems 2 and 3 in [ ]. Following a similar argument it should be possible to 
show that for sufficiently small fik the integral (88) is close to its approximation 
along the orbit of the Hill's problem. Such an argument is by no means trivial. 
A rigorous proof along these lines is in preparation and will be a subject of 
forthcoming work. For now we will quickly sketch how such an approximation 
works, but what follows should be regarded as a heuristic numerical support for 
the fact that for smaller Lik the integral (88) will be nonzero. This by no means 
is a proof of this fact. 

Applying the change of coordinates (10) to the integral (88) and passing to 
the limit with ll to zero, by direct computation (using (89) and (90)) we obtain 



k 




2 

3 

4 

5 

G 

7 

8 

9 

10 

11 

12 

13 



1.71124 
1.60503 
1.61216 
1.59042 
1.60485 
1.58405 
1.56759 
1.58011 
-1.5715 
1.58531 
1.56223 
1.5695 




(91) 
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where g Hlll (i) is the orbit on the one dimensional unstable manifold of the Hill's 
problem (11) at Lf m = (3- 1 / 3 , 0, 0, 3^/ 3 ), and where G™ 1 is given by 

with r 2 — x 2 + y 2 . 

To compute the integral M t HlU we need to specify the initial condition for the 
orbit q mn (t). This should correspond to the initial condition for the homoclinic 
orbit q°(t), which is at the intersection with the section {y = 0} (see Figure 
6 A.). Clearly q Hlll (t) will never cross such a section. Following [16], the initial 
conditions should lie on the tips of the "waves" of the orbit q mll {t). When con- 
sidering small fik for even k (recall that k is the number of waves which the orbit 
q°(t) makes) we should choose tips lying close to the zero velocity curve of the 
Hill's problem i.e. <? Hlll (0) = pk (see Figures 2 and 6 B.). For ^ with odd k we 
should choose tips g Hlll (0) = pt which lie far from the zero velocity curve. By 
choosing sufficiently large k the integral (91) provides an accurate approximation 
of (88) with sufficiently small /i^. The numerical results for this integral are 

ji /-Hill _ / -1.5464 for even k, (Q0 ] 

M * ~ \ 1.5464 for odd k, iy > 

which numerically supports the claim that for sufficiently small ^fe, for to = 0, 
we have a nontrivial zero of the Mclnikov integral (76) for the PRC3BP (See 
also Table 2). 

Remark 14- As mentioned before, a rigorous argument which follows the above 
outline is under preparation and will be the subject of forthcoming publication. 
Precise estimates which would prove the fact that §gf (0) tends to A/ t HlU are 
far from trivial. In the forthcoming work we also intend to provide a rigorous 
estimate for M t HlU instead of its numerical approximation (92). This will involve 
rigorous, computer assisted integration of (91), and is once again far from trivial. 



8. Transition Chains and Main Result 

In this section we will show that there exists a sequence of Lapunov orbits /(r^) 
for i = 1, . . . , N which survive for a sufficiently small perturbation e, and such 
that their stable and unstable manifolds intersect transversally 

W7. (n) rh WJ2 (r4+l) and W£ (r0 rhWi'. (n+l) . 

Such sequences are referred to as transition chains. Existence of such chains en- 
sures that we have an orbit which shadows the heteroclinic connections between 
the surviving tori generating rich symbolic dynamics (see [7] and [8]). Apart 
from this, from our argument it will also follow that for any surviving orbit l(r) 
we also have transversal intersection of its stable and unstable manifold 

W? Hr) rh Wf Hr) . 

This ensures that the chaotic dynamics of the PRC3BP which is implied by the 
existence of a transversal homoclinic orbit to l(r) (see Theorem 3 and Remark 
2) survives. 

We are now ready to rigorously reformulate our main theorem (Theorem 1). 
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B. y 




Fig. 6. The homoclinic orbit q° (t) in the (x,y) coordinates of the PRC3BP (A.) and in the 
Hill's coordinates (x,y) (B.). 

Theorem 11. For any fi from the sequence of masses {^tfc}fcL 2 f rom Theorem 2 
there exists a radius R(^) such that for sufficiently small eccentricities e of the 
elliptic problem there exists a Cantor set C C [0, R(fJ-)] such that for all r G <L 
the Lapunov orbits l(r) survive and are perturbed to invariant tori l e (r). What 
is more, there exists a transition chain l e (ri) for a sequence of radii < r± < 
r2 < ■ ■ ■ < r^r < R([i) such that r^ — ?T = 0(e 1 ^ 2 ). 

Remark 15. Let us note now that the proof of the above theorem will be based 
on some numerical results. First will be associated with the application of the 
KAM Theorem. For this we require that we have a twist property on the family 
of Lapunov orbits. This has been rigorously proved for sufficiently small ^ in 
Theorem 9, but the fact that we have twist for /Zfc with k = 1, 2, . . . has only been 
demonstrated numerically (see Section 4, Table 1 and Remark 9). Secondly, the 
argument requires that we verify that the Mclnikov function (76) has nontrivial 
zeros. This has only been demonstrated numerically in Section 7. 

We believe that the above can be verified using rigorous-computer-assisted 
methods. This is currently a subject of ongoing work. 

Proof (of Theorem 11). Let us fix a fi = By Theorem 9 and Remark 9 we 
have a radius R(fJ.) and a Cantor £ c [0, R{^)} of radii for which the Lapunov 
orbits survive the perturbation. In Section 7 we have verified that the Melnikov 
function (76) has nontrivial zeros. From Theorem 10 Remark 11 and Corollary 
1 it follows that by sufficiently reducing i?(/-i), we know that there exists a k > 
such that for radii r 1; r 2 G £ such that |ri — r 2 | < ne, for sufficiently small e we 
have 

W£(n) * W?(r 4 ) for ie{l,2}. 

We now need to show that we can find a sequence r\ < r 2 < . . . < rjv such 
that Ti € £ and rjy — ?T = 0{e 1 ^ 2 ), for which the gaps between r, and rj+i are 
smaller than ne. The existence of such a sequence follows from Remark 10. 
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9. Concluding remarks, future work 

In this paper we have shown that the chaotic dynamics observed for the planar 
restricted circular three body problem survives the perturbation into the planar 
restricted elliptic three body problem, when its eccentricity is sufficiently small. 
We have also shown that this dynamics is extended to include diffusion in energy. 
The diffusion proved in this paper covers a small range of energies. This is due 
to the fact that in our argument we use a Melnikov type method which does not 
allow us to jump between the "large gaps" between the KAM tori. An interesting 
problem which could be addressed is whether these large gaps can be overcome 
(this potentially could be done using techniques similar to [6] or [9]). 

Our result holds only for a specific family of masses of the primaries. The 
choice of these masses is such that they ensure the existence of the homoclinic 
orbit to L2, which is then us for the Melnikov argument. An interesting question 
is whether one can observe similar dynamics in real life setting, say in the Jupiter- 
Sun system. In such a case we will no longer have a homoclinic connection for 
the point L2. For the (circular) Jupiter-Sun system though we know that we 
have a transversal homoclinic connection for Lapunov orbits (see [13] and [23]). 
Such orbits could possibly be used for a similar construction. Our argument also 
required that we have sufficiently small eccentricities. It would be interesting 
to find out if the dynamics persists for the actual eccentricity of the Jupiter- 
Sun system. For this problem it is quite likely that applying the mechanism 
discussed in this paper would be very hard. Our argument relies on the use of 
the KAM theorem, which works for sufficiently small perturbations. To apply 
it for an explicit eccentricity seems a difficult task. Other methods could be 
exploited though. Instead of proving the persistence of the tori and trying to 
detect intersections of their invariant manifolds, one could focus on detection of 
symbolic dynamics for the diffusing orbits in the spirit of [23]. This seems a far 
more realistic target for the near future and is being currently considered as an 
extension of this work. 
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